Universality in edge-source diffusion dynamics 
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We show that in edge-source diffusion dynamics the integrated concentration N(t) has a universal 
dependence with a characteristic time-scale r = (A/V) 2 tt/ (AD), where D is the diffusion constant 
while A and V are the cross-sectional area and perimeter of the domain, respectively. For the short- 
time dynamics we find a universal square-root asymptotic dependence N(t) — N ^/t/r while in the 
long-time dynamics N(t) saturates exponentially at N . The exponential saturation is a general 
feature while the associated coefficients are weakly geometry dependent. 
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Concepts like diffusion and Brownian motion are cen- 
tral in a wide range of complex dynamical phenom- 
ena Q, Q including diffusion of ions through biological 
membranes, neutron diffusion in nuclear reactors, charge- 
carrier diffusion in semiconductors, diffusion of heat in 
any substance, diffusion of momentum in fluids, and dif- 
fusion of photons in the interior of the Sun 

0SHH In 

its simplest form of a scalar quantity c the diffusion dy- 
namics is governed by the linear partial differential equa- 
tion 



DV 2 c(r,i) = d t c(r,t), 



(1) 



where D is the diffusion constant. With the above nota- 
tion we have emphasized diffusion of matter (with con- 
centration c), but we note that the same equation also 
governs diffusion of energy such as in thermal problems 
where Eq. Q is often referred to as the heat equation. 
Despite its apparent simplicity the link between space 
and time variables r and t typically makes diffusion dy- 
namics strongly dependent on the geometry and initial 
conditions. In this paper we find an exception to this and 
report a universality in edge-source diffusion dynamics. 
Assuming perfect translation invariance along the z di- 




FIG. 1: A sketch of the cross section Q (white) with boundary 
(9fl (black curve) in the xy-p&ne of a rod with translational 
invariance in the z direction. At time t — 0, where the con- 
centration c(r, t) inside Q is zero, diffusion is suddenly turned 
on from the outside, where the concentration is kept at the 
constant value c (gray). 



rection, we consider a cross section f2 in the xy plane, see 
Fig. ^ where chemical species or heat is supplied at the 
boundary <9f2. For simplicity we imagine a situation with 
a constant concentration c outside Q while the domain 
itself is empty, c(r, t < 0) = for r € f2, before the onset 
of diffusion at t — 0. 

Recently, pressure-driven flow in steady-state was ana- 
lyzed in the framework of simple geometrical measures Q 
such as the cross-sectional area A = J n dr, the perime- 
ter of the boundary V — J gQ dl, and the compactness 
C = V 2 1 A. In the following we study the importance of 
these parameters for diffusion dynamics and in particular 
for the integrated concentration N(t) given by 



N(t) 



/ drc(r,<), 
Jn 



(2) 



with the limits N(t -> 0) =0 and N(t -> oo) = c A 
implied by the boundary and initial conditions for c. 

Dimensional analysis. As a first step in solving the 
dynamics of the integrated concentration N(t) one might 
estimate the time-scale t for filling up the domain. Ob- 
viously, increasing the area A results in an increasing 
filling time, while increasing the perimeter V or the dif- 
fusion constant D results in a decreasing filling time. By 
dimensional analysis we thus arrive at r - {A/V) 2 /D 
which, as we shall see, is indeed a good estimate since 
detailed analysis yields 



A 

T = [ v 



AD' 



(3) 



Short-time diffusion dynamics. On a short time-scale 
the diffusion is perpendicular to the boundary and thus 
the problem is quasi-one dimensional. By a short time- 
scale we here mean t <g; Lg n /D where L an is a character- 
istic length scale (such as the local radius of curvature) 
for shape variations along the boundary dfl. We base 
our analysis on the well-known method of combination 
of variables originally introduced by Boltzmann 0, 0]: 
normalization of length scales by V 'ADt reduces Eq. Q 
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TABLE I: The parameters a x Eq. © and ft Eq. JTTJ for the 
lowest eigenfunction for differently shaped cross sections. 



Cross section 



Circle 

Half-circle 

Quarter-circle 



1~ 



■nil 

(4+8/tt)^ 
2 

^2,1 
(4+16/ir) 2 



1.14 
- 1.08 
~ 1.00 



4/ 7 * ! ~ 0.69 
0.64° 
0.65° 



Ellipse(l:2) 


1.18' 




0.67" 


Ellipse(l:3) 


1.21' 




0.62" 


Ellipse(l:4) 


1.23' 


I 


0.58" 


Triangle(l:l:l) il 


tt 3 /36 ~ 


0.86 


6/tt 2 ~ 0.61 


Triangle(l:l:\/2)° 


5n 3 
16(2+v / 2) 2 


~ 0.83 


512/97T 4 ~ 0.58 


Square(l:l) 


tt 3 /32 ~ 


0.97 


64/tt 4 ~ 0.66 


Rectangle(l:2) 


5tt 3 /144 i 


rL 1.08 


64/tt 4 ~ 0.66 


Rectangle(l:3) 


5tt 3 /128 : 


± 1.21 


64/tt 4 ~ 0.66 


Rectangle(l:4) 


17tt 3 /400 : 


- 1.32 


64/tt 4 ~ 0.66 


Rectangle(l:oo) 


~ 7T 3 /16 1 


- 1.94 


64/tt 4 ~ 0.66 


Rectangle(w:h) 


7T 3 C-8 p 


4(fc + u>) 2 


64/tt 4 ~ 0.66 


Pentagon 


1.02' 




0.67" 


Hexagon 


1.05' 




0.68" 



"Data obtain ed b y finite-element simulations llll . 
See e.g. Ref. 111! for the eigenfunctions and eigenspectrum. 
c See e.g. Ref.ll2l for the eigenfunctions and eigenspectrum. 



to a ID ordinary differential equation with the solution 



erfc 



(4) 



where d± (r) is the normal-distance to the boundary and 
erfc(x) is the complementary error function. For t <C t 
the integrated concentration N(t) becomes 



N(t) « V 



ds c erfc( 



(5) 



where we have introduced the characteristic time scale 
defined in Eq. and N = cqA. The square-root de- 
pendence is a universal property for geometries with suffi- 
ciently smooth boundaries and dynamics deviating from 
this dependence is referred to as anomalous. Previous 
work on the heat content in the crushed ice model has 
reached results equivalent to Eq. © for the heat con- 
tents 0. 

Long-time diffusion dynamics. Since N(t) < N the 
result in Eq. JSJl is of course only meaningful for t <C 
r. When time becomes comparable to r a saturation 
will occur due to decreasing gradients in density. For 
structures without high symmetries the saturation will 
be accompanied by an onset of diffusion parallel to the 
boundary dVL, and in this limit the dynamics will be slow 
compared to the initial behavior, Eq. (J5J. To study this 
we first derive a continuity equation by applying Green's 



theorem to Eq. £[J, 

D [ d£n- Vc(r,£) = d t N(t). 
Jon 



(6) 



Here, n is a normal vector to dfl, the integral is a line 
integral along dQ, and — DVc is naturally interpreted 
as a current density. Next, we note that for long time- 
scales we have to a good approximation that n- Vc(r, t) cx 
c — N(t)/A is constant along the boundary so that 



N Q - N(t) cx d t N(t), 



(J) 



resulting in an exponentially decaying difference. This 
may also be derived from an eigenfunction expansion, 



c(r,t) 



= 1 



(8) 



which upon substitution into the diffusion equation yields 
a Hclmholz eigenvalue problem for <fi n and a n , 



(9) 



with <f> n (r) = for r 6 dil. Eq. JHJ) and the initial condi- 
tion c(r, t = 0) = imply that /n^ri( r ) = 1 an d thus 
fn = /n dr< M r )// n dr|<?!>„(r)| . The long-time dynam- 
ics is governed by the lowest eigenvalue yielding 



N, 



NCt) v-^ / t\ 

n 

1 -/^exp (- 4>r, (10) 



where 



fir, 



A cti 
A ' 



A 



off 



\S n dv<j) n {Y)\ 

Jn rfr K(r)| 



(11) 



As often done in optics , A^ can be interpreted as the 
effective area covered by the nth eigenfunction. Values 
for a selection of geometries are tabulated in Table 1. 
The circle is the most compact shape and consequently 
it has the largest value for (3 1 , or put differently the mode 
has the relatively largest spatial occupation of the total 
area. The normalized eigenvalue a Y is of the order unity 
for compact shapes and in general it tends to increase 
slightly with increasing surface to area ratio V j A. The 
modest variation in both a x and f3 1 among the various 
geometries suggests that the overall dynamics of N(t) 
will appear almost universal and that, e.g., Eq. ljlTJ|) for 
the circle (mathematical details follow below), 



N(t) 



4 

exp 

7oi 



( 2 f ( \ 



2 



5.783, (12) 



will account quantitatively well even for highly non- 
circular cross sections. 
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FIG. 2: A lin-log plot of N(t)/N Q as a function of t/r for 
a circular cross section (white circles) and for an arbitrarily 
shaped cross section (gray stars); both shapes are shown in 
the insets. The solid line shows the exact result for the cir- 
cle, i.e., the first 1000 terms in the infinite series Eq. (1141 . 
the dashed line shows the short-time asymptotic expression 
Eq. JSJj the dot-dashed line shows the long-time asymptotic 
expression Eq. Q12jl for the circle with a.-^ — j^/16 and 

Pi — 4/7oi> w hile the data-points are the results of time- 
dependent finite-element simulations. 



Analytical and numerical examples. In the following 
we consider a number of geometries and compare the 
above asymptotic expressions to analytical and numeri- 
cal results. For the numerics we employ time-dependent 
finite-element simulations ^(J an d solve Eq. (JTJ with a 
subsequent numerical evaluation of the integrated con- 
centration, Eq. J2J. 

The circular cross section serves as a reference and an 
illustrative example where we directly can compare the 
above limits to analytical and numerical results. Apply- 
ing the method of separation of variables yields || 

CcircO, t) = 1 _ 2 J ( 70: „r/a) e _^ nDt/a 2 ^ 



\ 7o,n J i(7o,J 

where a is the radius and r y m n is the nth zero of the mth 
Bessel function of the first kind, J m {lm,n) = 0- By a 
straightforward integration over the cross section we get 



N ciTC {t) 



1 



n=l '0,1 



exp 



2 

To,- 



7T 

16 



(14) 



where we have made the time-scale r explicit, see Eq. 
This result can also be derived using the continuity equa- 
tion, Eq. ©• For t > t the series converges rapidly, and 
keeping only the first term we arrive at Eq. Ijl2|l . 

In Fig. |21 we compare the asymptotic results, Eqs. (J5J 
and l|12|) with the exact result, Eq. (|14fl . as well as with 




FIG. 3: Plot of the deviation AN(t)/N from the circular 
case, Eq. I|15|l . as a function of t/r for different cross sec- 
tions (see insets). The data-points are results of finite-element 
simulations. The ellipses have eccentricities 1, 1.5,2 and 2.5, 
the circle sectors have angles 7r, 27r/3, 7r/2 and 7r/3, and the 
rectangles have aspect ratios 1,1.5,2 and 3. Note how the 
maximal deviation is less than ±0.03 or ±3%. 



time-dependent finite-element simulations. As seen both 
the short-time square-root and long-time exponential de- 
pendencies are in good agreement with the exact results 
as well as with the simulations. Fig. [3 also includes nu- 
merical results for an arbitrarily shaped cross section and, 
as suggested above, we see that Eqs. © and (tl"2l) account 
remarkably well even for this highly non-circular shape. 

In order to see how well Eq. Q14[l accounts for 
other non-circular geometries we have employed time- 
dependent finite-element simulations to numerically 
study the relative deviations AN/N from it, 



AN(t) _ N(t) - N ciic (t) 



(15) 



Figure |3 summarizes results for a number of geometries. 
In all cases the dynamics at small times t <C r is in full 
accordance with the predicted square-root dependence, 
Eq. JSJ , and for long times the predicted exponential de- 
pendence, Eq. I|12l) . fits the dynamics excellently. For 
the dynamics around t *~ r deviations from the circu- 
lar result are well within ±3% for the considered highly 
non-circular geometries. 

Discussion and conclusion. We have shown that edge- 
source diffusion dynamics in a rod of arbitrary cross sec- 
tion S7 has an intrinsic time-scale r = (A/V^tt/ (4D), 
with D being the diffusion constant while A and V are 
the cross-sectional area and perimeter of f2, respectively. 
Initially, the filling N(t) follows a universal square-root 
dependence N(t) — N y / T/r, irrespectively of the shape 
of the domain For longer times N(t) saturates expo- 
nentially at N . The saturation is governed by the low- 
est dimensionless eigenvalue a 1 of the Hclmholz equation 
rather than the full spectrum. Since a 1 depends only 
weakly on the geometry the dynamics becomes almost 
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FIG. 4: Plots as in Fig.|3]but here for two extreme geometries: 
a sector of a circular disc with angle 2 arctan(l/30) and a 
rectangle with aspect ratio 1/30, see insets. While these two 
shapes have nearly the same area and perimeter, the former 
fills slower than the circular cross section and the latter faster. 
The maximal deviations are —6% and +10%, respectively. 



universal. Numerically, we have observed that the devi- 
ation from strict universality is typically less than a few 
percent. 

The diffusion problem presented here relates to the 
question posed by Mark Kac 0] : " Can one hear the 
shape of a druml" . In the present diffusion problem 
knowledge about the short-time dynamics allows one to 
extract the area to perimeter ratio A/V while the shape 
itself cannot be inferred. For the long-time diffusion 
dynamics strict universality would require that different 
shapes have Helmholz eigenfunctions with the same set 
of eigenvalues {a n } (isospectrality) and effective areas 
{/3„}. However, since the answer to the question of Kac 
in most cases is positive, see however Refs. the 
eigenfunction properties {«„} and {P n } of different ge- 
ometries do differ. It is thus only the short-time dynamics 
which is strictly universal while, as mentioned above, the 
long-time dynamics depends weakly on shape through 
the first dimcnsionless eigenvalue a 1 and the correspond- 
ing dimcnsionless effective area (3 1 . 

Our simulations support these conclusions, see Fig. |3 
and even for extreme shapes such as the narrow disc sec- 
tor with angle 2 arctan(l/30) and the flat rectangle with 
aspect ratio 1/30 the deviations from Eq. (jT4f> are less 
than 10% around t ~ r, see Fig.0J These extreme shapes 
have almost the same area and perimeter while the eigen- 
function properties {«„} and {/?„} are very different, 
e.g., j3 l is constant 64/7T 4 for all rectangles, see Tabled 
while it scales as m -2 / 3 for the disc sector with angle 
7r/m. For any aspect ratio h/w the lowest eigenfunc- 
tion = sm(Trx/w) sin(7rt///i) of the rectangle is nearly 
uniformly distributed in Q, and the shape favors rapid 
perpendicular diffusion resulting in a filling slightly faster 
than for the circular shape. For the disc sector with angle 
ir/m the lowest eigenfunction <f> 1 = J m (7 m if /a) sin(m#) 



is confined to a region of width a/m 1 / 3 near the circu- 
lar edge. This shortens the effective perimeter resulting 
in a filling time longer than for the circular shape. It 
is thus possible to find extreme shapes where A/V is no 
longer the characteristic length-scale for the long-time 
diffusion dynamics. The sam e ap plies for a cross section 
with the shape of a tear drop 17] , and for more exotic ge- 
ometries, such as cross sections with a fractal polygonal 
boundary , we expect more severe deviations from the 
dynamics reported here. However, the deviations from 
strict universality obtained by extending the short-time 
scale r to the long-time regime are remarkably small. 

Apart from the fascinating and intriguing physics in- 
volved we believe our results are important to a num- 
ber of practical problems including mass diffusion in 
microfluidic channels and heat diffusion in arbitrarily 
shaped rods. 
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